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Multinucleon ejection model for Meson Exchange Current neutrino 
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A model is proposed to describe pairs or triples of nucleons ejected from nucleus as a result of Meson Exchange 
Current neutrino interaction. The model can be easily implemented in Monte Carlo neutrino event generators. 
It can provide a help in identifying true charge current quasielastic events and allow for better determination 
of the systematic error of neutrino energy reconstruction in neutrino oscillation experiments. 
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I. INTRODUCTION 

There is a lot of evidence for a significant multinucleon 
ejection contribution to the inclusive neutrino charge cur- 
rent (CC) cross section in the 1 GeV energy region 1 . 
On the experimental side, several recent nuclear target 
CCQE (CC quasi elastic) cross section measurements re- 
ported large values of the axial mass in disagreement 
both with older deuterium data and also with indepen- 
dent electroproduction arguments^. An explanation for 
the discrepancy would be that some events interpreted as 
CCQE are in fact due to very different dynamical mech- 
anism which typically leads to multinucleon emission. In 
fact, in the MiniBooNE (MB) experiment nucleons in the 
final state were not analyzed and only events with pions 
were identified and rejected. MB collaboration reported 
high statistics 2-dimensional muon inclusive cross section 
which is likely to include a large multinucleon ejection 
contribution and is a challenge for theoretical models^. 

On the theoretical side they are several computations 
which support the idea of the large multinucleon con- 
tribution. In the context of neutrino interactions the 
first was the Marteau model, based on the earlier ideas 
of Ericsson and Delorme^. The model uses the non- 
relativistic Fermi Gas in the Local Density Approxima- 
tion (LDA) approach. It includes elementary QE and 
A excitation interactions, Random Phase Approximation 
(RPA are in medium polarization effects important in the 
low Q 2 region^) and A width modification in the nuclear 
matter—. The model was later upgraded by Martini, Er- 
lcson, Chanfray and Marteau (MEChM model)!. With 
the inclusion of relativistic corrections MB 2D differential 
cross section can be reproduced^. 

Another approach to multinucleon ejection was pro- 
posed by Nieves, Ruiz Simo and Vicente Vacas 9 . The 
model is based on the earlier computations of the inclu- 
sive electron-nucleus cross section in the energy region of 
QE and A excitation peaks and also in the so called dip 
region between themi^. In the paper— the fit was done 
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to the MB 2D data with the best fit axial mass value 
equal to 1.08 ±0.03 GeV. 

Both above mentioned approaches are non-relativistic 
although several relativistic corrections were introduced 
based on the authors experience from dealing with the 
electron scattering data. 

An attempt to include relativistic effects in a more sys- 
tematic way was undertaken in a series of papers^. A 
2p-2h (two particles and two holes) contribution to both 
transverse and longitudinal responses coming from Me- 
son Exchange Current (MEC) and correlation diagrams 
were evaluated. Under some assumptions a relativistic 
MEC contribution to the neutrino inclusive cross section 
was also evaluated. 

The results of more rigorous computations done for 
electron scattering on lighter nuclei^ show that one can 
expect a large contribution to the transverse response 
coming from two-body current. But it seems to be im- 
portant to have a more realistic description of the nucleus 
ground state than it was assumed in all the previously 
mentioned papers. 

The situation is little confusing as the models which 
have better desription of the ground state, like the Spec- 
tral Function approach 14 , do not lead to better agreement 
with the MB data unless again a very large value of the 
axial mass is assumed^. It seems that it is important to 
have both Short Range Correlations (SRC) and MEC ef- 
fects in the same relativistic computational scheme which 
is however very difficult to achieve. 

Recently still one other effective approach to describe 
multinucleon contribution to the neutrino inclusive cross 
section was proposed. The Transverse Enhancement 
Model (TEM)i£ is based on the analysis of the electron- 
carbon scattering data and parameterizes the MEC effect 
as a modification of the magnetic electromagnetic form 
factor. The model predicts that the MEC contribution 
is less important at larger neutrino energies which can 
reconcile MB and NOMAD 17 axial mass measurements. 

It is clear that understanding of the multinucleon ejec- 
tion contribution is required for neutrino energy recon- 
struction which is in turn crucial for correct interpreta- 
tions of high precision neutrino oscillation measurements. 
Th^ftc?,!,^ contribution 
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which can be implemented in the neutrino Monte Carlo 
(MC) event generators 1 — and for that one needs also pre- 
dictions for the nucleons in the final state. With such a 
model one can start a direct comparison between what is 
measured and what is predicted in order at least of eval- 
uate the systematic error due to a lack of more precise 
multinucleon ejection model. Theoretical models provide 
only differential cross sections in the muon three momen- 
tum and a model of multinucleon ejection is something 
with which the theoretical models must be suplemented. 
Having in mind an importance of the predictions for the 
final state nucleons we propose a simple model to provide 
such information which can be confronted with the data. 



CCQE and np-nh on carbon with RPA 




neutrino energy [MeV] 



II. INCLUSIVE MUON CROSS SECTION MODELS 

The multinucleon ejection model of this paper is based 
on the information contained in the muon inclusive dif- 
ferential cross section only. If the kinematic variables are 
chosen to be energy and momentum transfer, one knows 
their joint probability distribution as a function of neu- 
trino energy. We will propose a procedure to predict 
nucleon momenta in the final state. 

We will discuss two models for which it is impossible 
to get exact predictions for the final hadronic system and 
one must rely on approximate modeling like the one we 
are discussing in this paper. 



A. Simplified Marteau (SM) model 

The first model we consider is a simplified version of 
the Marteau model. It was developed about 10 years ago 
and it includes RPA, A width, but not LDA effects. Its 
basic ingredients are presented in the reference^. Re- 
cently, elementary 2p-2h responses were added using the 
procedures described inl. 

Fig. rjjshow the comparison between predictions for the 
total cross section of the SM model with its full version. 
Here the term CCQE means that only one nucleon is 
ejected from nucleus. The SM model predicts larger cross 
section for the CCQE and lower for np-nh. Fig. [5] shows 
that for larger neutrino energies a part of the difference 
for QE is due to the fact that the SM model is more 
relativistic than the MEChM model. This is because in 
the computations the relativistic Lindhard function was 
used. 

The SM model predicts (like MEChM) both two- and 
three- nucleon ejection. The three nucleon contribution 
comes from the pionless A decay i.e. from reactions 
NA NN and NNA NNN with N standing for 
a nucleoid. 



FIG. 1. A comparison of predictions of the simplified Marteau 
(SM) model used in this paper and the full model called here 
MEChM. The target is carbon and all the nuclear effects 
are included. We show predictions for QE events (one nu- 
cleon ejection) and multinucleon (two or three) ejection. For 
MEChM the points are taken from 7 -. 



CCQE cross section on carbon; no RPA 



6 













SM model 






relativistic FG model 




*'.• 


MEChM model 



200 400 600 800 1000 1200 
neutrino energy [MeV] 

FIG. 2. A comparison of predictions for QE events before 
RPA nuclear effects. We show predictions of the SM model, 
MEChM model and also relativistic FG model. For MEChM 
the points are taken from 7 . 

B. Trensverse Enhancement Model (TEM) 

Ir>i& a new approach to describe CCQE scattering on 
nuclear targets is proposed. The model is easy to imple- 
ment in MC event generators. It is sufficient to modify 
vector magnetic form factors keeping all other ingredients 
of the CCQE model as in the free nucleon target case. 

The authors ofi£ proposed a universal transverse en- 
hancement function of Q 2 for the carbon target. For low 
Q 2 its form is determined by the scaling arguments while 
for high Q 2 (> 0.5 GeV 2 ) it is obtained as a fit to the 
inclusive electron cross section data from the JUPITER 
experiment. The prescription to include transverse en- 
hancement contribution in the numerical computations 
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np-nh cross section on carbon 
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FIG. 3. A comparison of predictions for the total multinucleon 
ejection cross section from three models: SM model, TEM and 
MEChM model. The target is carbon. 

amounts to the replacement: 

G p A ?(Q 2 ) GTO 2 ) = ^l + AQ^M-^G^iQ 2 ) 

(1) 

where G p ^ l (Q 2 ) are electromagnetic form-factors, A = 6 
GeV~ 2 and B = 0.34 GeV 2 . 

The most interesting feature of the TEM model is that 
it offers a possible explanation to the apparent contra- 
diction between low (MB) and high (NOMAD) neutrino 
energy Ma measurements: for energies up to ~ 700 MeV 
TEM predicts the CCQE cross section to be similar to 
CCQE with Ma = 1.3 GeV. For higher neutrino ener- 
gies the TEM cross section becomes significantly smaller 
and at E v ~ 5 GeV it corresponds to CCQE with 
M A ~ 1.15 GeV. 

As the TEM prediction for the MEC contribution we 
simply take the difference between the cross sections cal- 
culated with modified and standard magnetic form fac- 
tors: 



da/dra for np-nh on carbon, E v =700 MeV 
SM model — 
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Energy transfer [MeV] 

FIG. 4. A comparison of predictions for the differential cross 
section in energy transfer for multinucleon ejection from three 
models: SM model, TEM and MEChM model. In the case 
of MEChM model the results from old and new versions are 
shown separately. 
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FIG. 5. As in the previous figure, with different neutrino 
energy. For 1 GeV the predictions from old MEChM model 
are not available. 



d 2 a TEM d 2 a CCQE 



dqduj 



dqduj 



d 2 a cc Q E 
dqdu 



(G P M n ). 



(2) 



In the case of TEM model we assume that only two 
nucleon ejection takes place. 



C. Model comparison 

Fig. [3] shows the contributions to the overall cross 
section on the carbon target coming from multinucleus 
ejection. There is a significant difference in the predicted 



size of the effect. It is an important feature of the TEM 
model that for larger neutrino energies the MEC contri- 
bution becomes relatively small. 

Figs H] and [5] show differential cross section for the 
multinucleon ejection contribution as a function of energy 
transfer for two values of the neutrino energy: 700 MeV 
and 1 GeV. We show also the available results from the 
MEChM model. Note, that there are two versions of 
the MEChM model with different elementary 2p-2h re- 
sponse functions. We followed the procedures described 
in 7 and obtained new responses as functions of the vari- 
able x = Q 2 /2Muj — (q 2 ~ lo 2 )/2Mlo where uj and q 
denote energy and momentum transfers and M is nu- 
cleon's mass. However, the new response functions lead 
to rapidly increasing multinucleon ejection cross section 
as a function of neutrino energy. The behavior seems to 
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be unphysical and in this paper we rely on the original 
Marteau responses. In Figs [U and [3] one can see signifi- 
cant differences between the models which will translate 
into predictions for the multinucleus ejection. The old 
MEChM model predicts the pronounced peak at low en- 
ergy transfer which is probably non-physical and comes 
only from oversimplified assumptions for the response 
function. The SM model predicts the peak which is 
smaller in size and for larger values of the energy transfer 
develops a shape which is close to the new version of the 
MEChM model. 

In the case of TEM the sharp fall down of the differ- 
ential cross section at ~ 410 MeV comes from the way 
in which the model was implemented. We followed the 
original paper and assumed the target nucleon to be at 
rest. Pauli blocking effects are introduced by means of 
the Q 2 depending suppression function. 



III. NUCLEON EJECTION MODEL 

Two basic assumptions of the model is that the en- 
ergy and momentum are transfered to two (or three) nu- 
cleons simultaneously and that there are no correlations 
between inital state nucleons. The second assumption 
can appear to be an unrealistic but its impact on the 
final state predictions is rather mild. 

The scheme of the procedure to generate nucleon final 
state is as follows: 

• two (or three) nucleons are selected from the Fermi 
sphere of radius 220 MeV (it is assumed that inter- 
action occurs on carbon) 

• four momentum of the hadronic system is calcu- 
lated by adding four momenta of selected nucleons 
and energy and momentum transfered by the inter- 
acting neutrino 

• Lorentz boost to the hadronic center of mass sys- 
tem is done 

• two (or three) nucleons are selected isotropically in 
the hadronic center of mass system 

• boost back to the laboratory frame is performed. 

The energy balance is done based on the assumptions 
that inital state nucleons are in the potential well of the 
depth V = Ef + 8 MeV (Ef is the Fermi energy). 

• Fermi energy is subtracted from each initial state 
nucleon 

• for each nucleon in the final state (in the LAB 
frame) the energy is reduced by the amount of 
8 MeV adjusting its momentum so that they re- 
mains on-shell. 



The above procedure allows for a smooth distribution 
of nucleon momenta in the final state. 

Pauli blocking was not imposed for two reasons. First 
of all, the final nucleon state is calculated for events about 
which it is known that they occur and a reduction of the 
cross section is not required. Secondly, a smooth spec- 
trum of final state nucleons should be obtained. Clearly, 
other prescriptions to account for the energy balance can 
be invented but the above procedure seems to the sim- 
plest one and satisfying basic physical requirements. 

It is an interesting observation that the above de- 
scribed algorithm actually introduces some correlations 
between initial state nucleons. Not all initial state nu- 
cleon configurations give rise a hadronic system which in 
the center of mass frame has the invariant mass larger 
than 2M (or 3M for three nucleon ejection; M is nu- 
cleon's mass) . In the numerical code initial state nucleon 
configurations are selected until acceptable nucleon mo- 
menta are found. 

We investigated an impact on the model predictions 
if a pair of nucleons in the initial state is chosen to be 
always in the deuteron-likc configuration i.e. with oppo- 
site three-momenta. It turned out that the results arc 
virtually unchanged. The most important assumption 
of the model is that of the isotropic distribution of final 
nucleons in the hadronic center of mass frame. 

Clearly, in the MC event generator there should also 
be included a contribution from SRC large momentum 
nucleon pairs. This contribution can be obtained in the 
SF formalism with a straigthforward MC implementation 
(a simpler alternative is the Bodek- Ritchie model^). In 
the SF approach the lepton interacts with only one nu- 
cleon which absorbs all the energy and momentum that 
is transfered to the hadronic system. The other nucleon, 
with approximately opposite initial three momentum is 
a spectator only. 

In the MC simulation all the above described steps 
should be followed by Final State Interactions (FSI). 



IV. RESULTS 

The series of figures show the comparison of predictions 
from both models for the inclusive muon cross section 
(TEM and SM model) with the same model of multin- 
ucleon emission. All the distributions are normalized to 
the same area so that only their shapes are compared. 



A. Total kinetic energy 

Figs [5] and [7] show the total kinetic energy of all the 
nucleons in the final state. Even if individual nucleons 
cannot have reconstructed tracks due to too small mo- 
menta, the total kinetic energy can be measured as a 
vertex activity. One can see that the predictions from 
two models differ significantly. The SM model predicts 
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FIG. 6. Total kinetic energy of the final state nucleons. Neu- 
trino energy is 800 MeV. 
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FIG. 7. Total kinetic energy of the final state nucleons. Neu- 
trino energy is 1500 MeV. 



more events with large energy deposit. In order to un- 
derstand this one should look back at the Figs [4] and [5] 
TEM predicts much smaller cross section at large values 
of the energy transfer which translates to the total kinetic 
energy of the ejected nucleons. 

In the case of the SM model we would like also to show 
separate contributions from two- and three- nucleon final 
states, see Fig. [5J It can be seen that three nucleon ejec- 
tion contribute significantly to the maximum at larger 
values of the kinetic energy. 



B. Nucleon momenta 



FIG. 8. Contributions from two- and three- nucleon final 
states to the overall SM model predictions for the total kinetic 
energy of ejected nucleons. Neutrino energy is 1500 MeV. 
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FIG. 9. Momentum of the most energetic nucleon in the final 
state. Neutrino energy is 800 MeV. 



and SM model is that the first one predicts a maximum at 
the value ~ 500 MeV/c while the other model prediction 
have the peak cut off and the extended large momentum 
tail. 

Figs[TT1and[T2lshow the distributions of the momentum 
of the second most energetic nucleon in the final state. 
In the case of two nucleon ejection it is also the lowest 
momentum nucleon but in the case of the SM model there 
is also a fraction of events with three ejected nucleons. 
One can be surprised that the second nucleon can be also 
quite energetic. 



Figs IH1 and [TU] show the distributions of the momenta 
of most energetic nucleon in the final state. This is an 
important prediction because it can tell us how likely 
it is that a MEC event will have a nucleon above the 
detection threshold. The main difference between TEM 



C. Angular distribution 

We calculated the cosine of the angle Oqe between the 
most energetic nucleon momentum and the hypothetical 
three momentum of the nucleon if the interaction were 
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FIG. 10. Momentum of the most energetic nucleon in the 
final state. Neutrino energy is 800 MeV. 



Second nucleon momentum (E v =800 MeV) 





0.12 - 




0.1 






E 


0.08 - 




>, 




co 


0.06 


5 




co 


0.04 




0.02 - 











SM model 












T! 


EM - 





















































































100 200 300 400 500 600 700 800 900 
momentum [MeV/c] 

FIG. 11. Momentum of the second most energetic nucleon in 
the final state. Neutrino energy is 800 MeV. 
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FIG. 12. Momentum of the second most energetic nucleon in 
the final state. Neutrino energy is 1500 MeV. 



FIG. 13. Cosine of the angle between the most energetic 
nucleon three momentum and hypothetical three momentum 
of the nucleon if the interaction were CCQE and the target 
nucleon at rest. Neutrino energy is 800 MeV. 
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FIG. 14. Reconstructed neutrino energies for MEC interac- 
tions of E v = 1500 MeV neutrino. Contributions from two 
and three nucleon ejection are also shown separately. 



CCQE with the target nucleon at rest {9qe is calculated 
based only on muon's kinetic energy and production an- 
gle). In experimental analyses 8qe is often used to se- 
lect CCQE enriched samples of events with a criterium 
defined as Qqe < 30°. Fig. [13] show that for neutrino en- 
ergy 800 MeV the distribution is strongly forward peaked 
and the nucleons ejected due to MEC mechanism seem 
to mimic those which arise after genuine CCQE interac- 
tion. For neutrino energy 1500 MeV the results are very 
similar. 



D. Energy reconstruction 

Fig. [HI shows the distribution of E®J^ obtained within 
the SM model. The true neutrino energy is 1500 MeV. 
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FIG. 15. ELAMM predictions for the total multinucleon ejec- 
tion cross section is compared with: SM model, TEM and the 
MEChM model. 



FIG. 16. ELAMM predictions for the total kinetic energy 
of the final state nucleons is compared with: TEM and SM 
model. Neutrino energy is 800 MeV. 



is defined in the standard way by assuming that 
only final muon is detected, the interaction was CCQE 
and the target nucleon was at rest. Two peaks seen on 
the figure come from various ingredients of the SM model. 
The most important observation is that on average the 
values of E® e f are much smaller than the true neutrino 
energy. 



V. DISCUSSION 
A. ELAMM 

In order to get still better intuition how large uncer- 
tainties of the model can be, we would like to introduce 
still another effective MEC 'model', the one defined by 
the large axial mass itself! In the 'model' which we shall 
call ELAMM (Effective Large Axial Mass Model): 
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FIG. 17. ELAMM predictions for the total kinetic energy 
of the final state nucleons is compared with: TEM and SM 
model. Neutrino energy is 1500 MeV. 



d 2 a MEC d 2 a CCQE 



dqduj 



dqduj 



(M A = 1.35) 



B. Isospin 



d 2 a CCQE 

dqduj 



(M A = 1-05) 



(3) 



ELAMM has an important advantage that it repro- 
duces very well the MB 2D differential cross section data. 

Figs mi mi and LIZ] show that ELAMM predictions are 
quite different from the TEM and the SM model. It can 
be surprising that at the typical MB neutrino flux ener- 
gies ~ 700 MeV the ELAMM multinucleus ejection cross 
section is much smaller than that of the SM model. Also 
ELAMM np-nh cross section increases for larger neutrino 
energies which is not the case for other models. It is 
rather clear that the way in which the ELAMM describes 
the excess of events due to MEC dynamics is rather dif- 
ferent to the one provided by TEM and SM model. 



The question which was avoided so far but which is 
very important from the experimental point of view (it 
is much easier to detect protons than neutrons) is that 
of isospin state of ejected nucleons. 

In the case of SRC nucleon pairs it is known that it 
is much more likely to have proton-neutron pairs than 
neutron-neutron or proton-proton ones^j22. It is not ob- 
vious in which way these results should be included in the 
model predictions. All the available microscopic models 
of MEC depart from the (local) Fermi Gas model. Purely 
combinatoric arguments (number of possible p-n and n-n 
pairs) lead to different conclusions than SRC arguments 
(p-n pair is six or more times more likely than n-n) . In the 
case of SM model (as it is also for the full MEChM model) 
an extra complication comes from the three nucleons fi- 
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FIG. 18. Distribution of total ejected nucleons kinetic energy 
averaged over MB flux. 



FIG. 19. Distribution of most energetic nucleon momentum 
averaged over MB flux. 



nal states for which it is even more difficult to make any 
assumptions. Altogether we think as reasonable to ex- 
pect that in CC neutrino reactions ~ 80% of final states 
nucleons are protons and only ~ 20% are neutrons. For 
CC antineutrino reactions the isospin composition of the 
final state is the opposite with ~ 80% of neutrons and 
~ 20% of protons. 

TEM gives us no hints about final nucleons isospin 
states and we would like to propose the same rules as 
those assumed for the microscopic SM model. 



C. Flux integrated predictions 

We would like also to show some flux averaged predic- 
tions. In the case of SM model results should be treated 
with a caution because the model is basically a non- 
relativistic one (even if some relativistic corrections are 
included) and neutrino fluxes in experiments like MB, 
T2K, MINERvA extend to the energies at which rela- 
tivistic effects are important. The minimal consistency 
requirement that we checked is that at large energies be- 
haviour of the SM is stable. 

Figs HH] and 033 show predictions from three models 
(including ELAMM) for the total kinetic and maximal 
nucleon momentum integrated over the MB flux. As 
expected, the differences are quite large and they can be 
understood by the previous disussion. 



VI. CONCLUSIONS 

In the paper a multinucleus ejection model is proposed 
which can be easily implemented in neutrino MC event 
generators and can help in selecting true CCQE events. 
The model proposed in the Section IIIII can be used in 
conjunction with any model providing predictions for 
the charge lepton inclusive neutrino-nucleus cross sec- 



tion. The models of muon inclusive cross sections dis- 
cussed in the paper are only approximations and it can 
be shown that they do not reproduce MB double differ- 
ential cross section too well2!; the problem is that there 
are no other available (in the form applicable to MC sim- 
ulations) models to predict muon inclusive cross section 
due to MEC mechanism. Our model is very simple but in 
reality, there are also FSI effects which must be included 
and which would smear out many details of a better jus- 
tified multinucleon emission model. 

The quantitative understanding of the size and kine- 
matical characteristics of the multinucleon ejection con- 
tribution can only come from coincidence charged lepton 
and final hadronic system measurements in several ex- 
clusive channels. This is clearly a very demanding goal. 

The observable which in the minimal degree depends 
on the details of a multinucleon ejection model is the 
integrated kinetic energy deposit of nucleons. In order 
to compare with the model predictions for individual nu- 
cleon momenta a detector with a low threshold for recon- 
struction proton tracks is a preferable choice. The liquid 
argon target with the threshold as low as ~ 300 MeV/c 
seems to be the best one. 800 MeV neutrinos can produce 
a sizable fraction of MEC events with two reconstructed 
proton tracks. 
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